Direct measurements of hydrophobic slippage using double-focus fluorescence 

cross-correlation. 
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We report results of direct measurements of velocity profiles in a microchannel with hydrophobic 
and hydrophilic walls, using a new high precision method of double-focus spacial fluorescence cross- 
correlation under a confocal microscope. In the vicinity of both walls the measured velocity profiles 
do not turn to zero by giving a plateau of constant velocity. This apparent slip is proven to be 
due to a Taylor dispersion, an augmented by shear diffusion of nanotracers in the direction of flow. 
Comparing the velocity profiles near the hydrophobic and hydrophilic walls for various conditions 
shows that there is a true slip length due to hydrophobicity. This length, of the order of several 
tens of nanometers, is independent on electrolyte concentration and shear rate. 
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PACS numbers: 82.70.Dd, 83.80.Qr, 82.70.-y 

For more than hundred years scientists and engineers 
have assumed and successfully applied no-slip boundary 
conditions to model experiments in fluid mechanics 
However, it has been recently well recognized, that the 
success of this famous no-slip postulate reflected mostly 
a macroscopic character and insensitivity of old exper- 
iments. Reducing the size of investigated systems to 
micro- and, especially, nanodimension led to a very defi- 
nite conclusion that the no-slip condition does not always 
apply 0- It is now clear that many systems should al- 
low for an amount of slippage, described in terms of a 
slip length: v s = bd z v, where v s is the slip (tangential) 
velocity at the wall and the axis z is normal to the sur- 
face. What, however, remains a matter of active debates 
is the amplitude of slip, and its variation with interfacial 
properties and parameters of the flow. 

From the theoretical point of view the situation is 
reasonably clear. Slippage should not appear on a hy- 
drophilic surface, except as at very high shear rate [3(. 
A slip length of the order of hundred nanometers or 
smaller is, however, expected for a hydrophobic sur- 
face 0, [H, Q. On the experimental side, no consen- 
sus is achieved so far. While some experimental data 
are consistent with the theoretical expectations both 
for hydrophil ic B , H, & EH 11 1 and hydrophobic sur- 
faces 0, la HI E37 some other reports completely escape 
from this picture with both quantitative (slippage over 
hydrophilic surface, shear rate dependent slippage, rate 
threshold for slip, etc) and quantitative (slip length of 
several /ims) discrepancies (for a recent review see [HI). 
Clearly, in order to rationalize the experimental situa- 
tion, new data are necessary. These data should prefer- 
ably be obtained with a new experimental technique. 

Basically, two types of experimental methods were 



used to study boundary conditions. High-speed force 
measurements performed with the surface forces appa- 
ratus (SFA) [a III or atomic force microscope (AFM) 
allows to deduce a drag force, with the subsequent com- 
parison with a theory of a film drainage This ap- 
proach, being extremely accurate at the nanoscale, does 
not provide visualization of the flow profile, so that this 
type of measurements should be identified as indirect. 
Direct approaches to flow profiling, or a velocimetry, 
take advantage of various optics to monitor tracer parti- 
cles [l5, 1(1 17 1 . Their accuracy is normally much lower 
than that of force methods due to relatively low optical 
resolution, system noise due to polydispersity of tracers, 
and difficulties in decoupling of directed flow from diffu- 
sion. As a consequence, it is normally expected that a 
slippage of the order of a few tens nanometers cannot be 
detected by a velocimetry technique. 

In this Letter we report direct high-precision measure- 
ments at the nanoscale performed with a new optical 
techni que. As an alternative to the existing FTIR [HI], 
M-PIV and TIRV methods we here use a re- 

cently suggested technique, based on a double-focus spa- 
tial fluorescence cross-correlation (DF FCS) [Tij . Our 
method allows one to use much smaller tracers, an or- 
der of magnitude higher shear rates, and to get orders of 
magnitude better statistics as compared with the state of 
the art. These give us at least an order of magnitude im- 
provement in accuracy compared with other velocimetry 
approaches. Results obtained for various experimental 
conditions allow us to deduce the true hydrophobic slip 
length, which is proven to be of the order of several tens 
of nanometers, and is independent on electrolyte concen- 
tration and shear rate. 

Our microchanel was formed by a three-layer sand- 
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FIG. 1: Schematic of the experimental setup. Abbrevi- 
ations are: MC - MicroChannel, WP - Wollaston Prism, 
DBS - Dichroic Beamsplitter, EF - Emission Filter, APD - 
Avalanche Photodiode, and PH - Pinhole. 



wich construction. The lowest layer was a standard mi- 
croscope cover slide made of borosilicate glass with a 
thickness of 170 /im, a root-mean-square roughness of 
the range 1-2 nm. The water advancing contact angle 
on this slide was measured to be below 5°. The chan- 
nel itself was created by cutting out a hole in an adhe- 
sive polymer film (Tessa, Germany) with a thickness of 
around 100 /im, that forms the smallest dimension of the 
channel, directed along the z axis. The channel exten- 
sion along the y axis (along the wall and perpendicular 
to the flow direction) is about 1.5 mm. Finally, the top 
layer was formed by a 1-mm-thick cover glass. Its sur- 
face was made hydrophobic by silanization and the water 
contact angle was measured to be 85 — 90°. An opti- 
cally transparent polycarbonate block served as a sup- 
port and for connection of the chamber to the external 
flow system. A hydrostatic pressure gradient was created 
by a system of two beakers at different heights, which al- 
lowed us to vary a shear rate near the wall in the range 
A = 800 — 3000 s -1 . As a tracers we used fluorescently 
labeled latex spheres, carboxylate-modified FluoSpheres 
580/605 (Molecular Probes, Eugene, Oregon). The par- 
ticles had a radius of R ~ 20 nm and a polydispersity of 
about 20%. Experiments were carried out in water and 
NaCl aqueous solutions with concentrations in the range 
between 10~ 5 mol/L and 10~ 2 mol/L. 

The scheme of the DF FSC method is shown on FigQ] 
and is described in details before [HI . Briefly, we used a 
commercial FCS setup (Carl Zeiss Jena, Germany) con- 
sisting of the module ConfoCor 2 and the inverted micro- 
scope model Axiovert 200. For the present experiments, 
we employed a water immersion objective (Zeiss, C- 
Apochromat 40x , NA 1.2). The optical system was mod- 
ified so that an external laser beam could be coupled into 
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FIG. 2: Schematics of the basic idea of the double- focus spa- 
tial fluorescence cross-correlation method. Two laser foci are 
placed along the x axis separated by a distance of 6 /im. They 
independently record the time-resolved fluorescence intensi- 
ties Ii(t) and 12(f). The forward cross-correlation of these 
two signals yelds G(t). Two foci are scanned simultaneously 
along the z axis to probe the velocity profile v(z). 



the confocal optics. For fluorescence excitation, the 543- 
nm line of a 1-mW helium-neon laser was used. The laser 
beam was split by means of a Wollaston prism. Behind 
the prism, the two beams are polarized perpendicularly 
to each other and exhibit an angular separation of 0.5°. 
After passing through two additional lenses these beams 
are fed into the confocal microscope. Our alignments re- 
sult in two optically equivalent, almost diffraction- limited 
laser foci (diameter ~ 400 nm, height ~ 3 /im) separated 
by a distance of 6.0±0.1 ^m in object space as is schemat- 
ically shown in Fig. [2j As the fluorescence tracers are 
flowing along the channel they are crossing consecutively 
the two foci, producing two time-resolved fluorescence in- 
tensities I\ (t) and I2 (t) recorded independently from the 
avalanche photo-diodes APD 1 and APD 2. The time 
cross-correlation function of h(t) and hf) can be cal- 
culated as G(t) = {h(t)I 2 (t + T)) t / {h(t)) t (I 2 (t)) t and 
typically exhibit a local maximum. The position of this 
maximum tm is characteristic of the local velocity of the 
tracers. 

To determine the velocity profile we have scanned the 
foci position across the channel. At each z position, a se- 
ries of 10 independent data acquisitions was carried out. 
The acquisition time was either 30 s or 60 s, necessi- 
tating longer measurements close to the channel walls, 
where small flow velocities are found. Indeed, consider 
the worst case when the foci are centered on the wall. 
Since the concentration in particles is about 1 per fem- 
tolitcr, the number of particles carried by the shear flow 
which enter the focus half-elliptical area during 60 s is 
about N = 6 x 10 4 . This gives a satisfactory signal to 
noise ratio y/~N larger than 10 2 . The independent cross- 
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FIG. 3: Typical velocity profile u(z) measured in a ~ 100 fim 
channel with 10 -4 mol/L NaCl solution. Inset: The observed 
velocity profile in the vicinity of the wall and schematics of 
the procedure for a determination of the apparent slip length 
(frapp ~ 445 nm) and the wall location (z = 50.7/zm). 

correlation functions acquired at position z were fitted 
with a Gaussian function for more precise determination 
of tm(z), yielding the particle velocity v(z) = As/tm(z), 
where As is the distance between foci. For every salt con- 
centration we have repeated experiments several times 
with freshly prepared channels and varied the pressure 
gradient. 

A typical measured velocity profile is shown in Fig. [3] 
As expected, the central region, where the velocity of 
the tracer particles reflect that a liquid [3], the profile 
exhibits the parabolic shape predicted by the classical 
theory. However, when foci presumably enter the wall, 
this parabola does not turn to zero by giving a plateau 
of non-zero constant velocity. This is observed for both 
hydrophobic and hydrophilic surfaces. The apparent ve- 
locity at the plateau region is always higher for a hy- 
drophobic surface and decreases with added salt and as 
shown in Fig. |4j 

The apparent velocity at the wall is too large to reflect 
the true liquid slippage over it. Earlier estimates 
suggested that in the vicinity of the wall the tracers are 
submitted to a Taylor dispersion, e.g. their diffusion is 
augmented by shear, enhancing a migration speed in the 
direction of flow. Now we model this effect precisely. 
Like in [Til we assume an crgodic system and interpret 
the time cross-correlation G(r) as: 

= //»i(r^ 2 (rO$( r , r ' ;T )d 3 rdV (1) 
C 2 fji 1 (r)i 2 (r')d 3 rd 3 r' 

where r = (x, y, z), r' = (x 1 ', y', z') and the average con- 
centration of labelled particles is denoted by C. The 
real-space detection efficiencies ii(r) and i2(r) for fo- 
cus 1 and 2 were given as ellipsoidal Gaussian functions. 
The function <!> is given as the solution of the advection- 
diffusion equation from a point source in a flow field with 
uniform velocity. In the vicinity of the wall, the parti- 
cles are repelled by an electrostatic force, F, so that they 



FIG. 4: The apparent velocity at the plateau region, «a PP , 
at the hydrophobic wall (triangles) and hydrophilic wall (cir- 
cles) as a function of concentration of NaCl. Concentration 
10~ 6 mol/L corresponds to a case of pure water. Dashed and 
solid curves show the values predicted for b = 100 nm and 
6 = 0, correspondingly (A = 1750 s _1 ). Inset shows typical 
calculated isoconcentration lines {c(x, z) =Ct) distorted by 
dispersion together with isolines of light intensity profile at 
the downstream focus {iz(x, z) =Ct). 

do not fill up completely the ellipsoidal Gaussian light- 
ened region. Keeping the same notation for the light 
intensity, we thus replace i\{r) by i\(r)c e (r), where 
c e (r ) is the particle equilibrium concentration profile at 
the upstream focus, to be detailed below. Wc also take 
into account the velocity gradient in the advection flows, 
thereby introducing the mechanism for Taylor dispersion. 
The advection-diffusion equation to be solved is: 

dc T + (Xz + b)dc x + d(wc) z = DV 2 c (2) 

where w = F j (6irRfi) is the migration velocity of a par- 
ticle along z in a fluid of viscosity (i and D is the Einstein 
diffusion coefficient. Here we neglect the hydrodynamic 
interactions between the particles and wall. Since the dis- 
tance from the wall h 3> R, the energy of electrostatic in- 
teraction of a particle with the wall is U = q<pi exp(— kH), 
i.e. the particle is considered as a point charge q [Corre- 
spondingly, F = —dll/dh]. Here 4>i is the surface poten- 
tial of the wall at the given concentration of electrolyte 
with an inverse Debye length k. The charge is given by 
q = A-KR 2 q s = iirReoefoinR + 1), where q s is the surface 
charge density, and is the surface potential of the par- 
ticle. Then, according to the Boltzmann law, the equilib- 
rium concentration of particles in the vicinity of the wall 
is c e (r) = Co exp(— U/ksT) = cq exp(— A exp(— re/i)), 
where A = Air RfafatoeCL + nR)/(kBT). For </>i we used 
data of data for </>2 were smaller according to 

electrokinetic measurements. With these parameters the 
values of A were of the order of 25. Instead of solving 
for an initial point source, multiplying by i\{r )ci(r ) and 
calculating the integral on r like in the numerator of EqQJ 
we may from linearity of Eq[2] solve for the initial cloud of 
illuminated particles with concentration i\(r )c\(r ). The 
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FIG. 5: Concentration at upstream focus from balance of 
Brownian motion and electrostatic repulsive forces for the 
typical case of 10 -4 mol/L NaCl, and intensity at the up- 
stream focus when centered at 100 nm from the wall. The 
concentration of illuminated particles is the product of these 
quantities. 



result c(r, r) is then multiplied by i^ij ') and integrated. 
This integral goes through a maximum at some time tm , 
which is interpreted as the transit time between foci, like 
in Fig. [21 Typical distributions of the concentration c e 
and the intensity i\ are shown in Fig. [5l 

For simplicity of the numerical analysis, we have re- 
duced the 3D problem to a 2D one by integrating Eq. [2] 
along y. Then we obtain for c = cdy an equation 
of the same form, with now V 2 = 9 2 2 + d^ 2 ■ The initial 
condition also is similar. The non-penetration condition 
d z c = applies automatically at the wall as a result of 
the repulsive electrostatic force. The equation advection- 
diffusion with these initial and boundary conditions was 
solved with a commercial finite elements software (COM- 
SOL). Experimental values of the shear rate were used. 
The calculation was repeated for several values of b in the 
interval from zero to 100 nm and for several values of the 
distance of the foci to the wall. In each case, the calcu- 
lated tm provided the apparent transfer velocity between 
foci, As/tm- 

Our calculations allow one to reproduce both the ve- 
locity profile close to the wall and the plateau region (see 
the solid dots in the inset of Fig. [3]) . By superimposing 
these two parts of the velocity profile for a given concen- 
tration of salt, we unambiguously determine a position 
of the wall in the experiment. An apparent slip length 
was then obtained by fitting a straight line through the 
points of the velocity in term of the distance to the wall, 
similarly to the inset in Fig.[3l The apparent slip lengths 
is about 445 nm for all electrolyte solutions, but is much 
larger in pure water, being of the order of 2 /im. These 
results slightly differ from previously reported [lit due 
to a different, not supported by our current model, way 
of determination of the wall position. As become evident 
from Fig. [3] the fit of experimental data always gives the 
values of the apparent slip close to predicted by the the- 
ory. Thus, we conclude that the large observed value 
of the apparent slip at the hydrophilic wall arc fully at- 



tributed to a Taylor dispersion affected by electrostatic 
interaction of nanotracers with the wall. 

The velocity profiles calculated for a slip wall are con- 
sistent with data obtained near hydrophobic surfaces, 
and the apparent slip at the hydrophobic wall is found to 
be 60 — 70 nm larger than in case of a hydrophilic wall for 
all salt concentrations. It follows from our modelling that 
the contributions of the Taylor dispersion for the no-slip 
and slip are of the same order, so that the difference of the 
apparent slip lengths is close to the actual ones. Alterna- 
tively, the true hydrophobic slip length was determined 
by comparison the apparent velocities at plateau regions 
(see Fig. |4]). This procedure is more accurate since it does 
not suffer from possible errors introduced by fitting and 
even does not depend of the choice of the wall location. 
Fig. 2] shows that the results for hydrophobic surfaces are 
always below the curve computed for b = 100 nm. Note 
that the presence of electrolyte has little effect on the 
value of the hydrophobic slip length. Neither found we a 
dependence of measured values on the shear rate. In the 
numerical and experimental examples we used here as an 
illustration of our approach the shear rate close to the 
wall (1750 s _1 ) was larger than in [IF 



17J and compara- 



ble to [12j . There are indications that shear rate strongly 



influence the value of apparent slip in all range of shear 
rates we used, but the true hydrophobic slip however re- 
mains the same as in the discussed data. 



In summary, we have performed an experimental study 
of a flow of water-electrolyte solutions in microchanels by 
using a new velocimetry technique. Our experiment is 
in favor of no-sli p b oundary conditions for a hydrophilic 
surface @, H, foL llOj. It is very unlikely that there ex- 
ists some minimal slip over hydrophilic surfaces as sug- 
gested before [l2, 17 1. We have also demonstrated that 
there is no possibility that flow exhibits a hydropho- 
bic slip length larger than 80-100 nm. Therefore, the 
slip effect is not as extreme as many authors have re- 
ported [l3l . flij . However, it is quite large if we con- 
sider a simple molecular model for slip at our contact 
angle [H, which might be a good indication to a two- 
layer model @, (b/6 « n/iig « 50, where 5 is the 
thickness of the adjacent to the hydrophobic surface gas 
layer with viscosity fi g , which suggests S m 1 — 2 nm). 
Essentially, our DF FCS approach allowed us to use very 
small particles, to reach a very large shear rate, and to re- 
duce dramatically an error in measurements due to orders 
of magnitude better statistics than known methods 21 1 . 
We thus believe our Letter concludes the discussion about 
boundary conditions at hydrophilic and hydrophobic sur- 
faces. 
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